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Pumping in quantum dots and non-Abelian matrix Berry phases 
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We have investigated pumping in quantum dots from the perspective of non-Abelian (matrix) 
Berry phases by solving the time dependent Schrodinger equation exactly for adiabatic changes. 
Our results demonstrate that a pumped charge is related to the presence of a finite matrix Berry 
phase. When consecutive adiabatic cycles are performed the pumped charge of each cycle is different 
from the previous ones. 
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I. INTRODUCTION 
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FIG. 1: Adiabatic path in the parameter space is shown. An 
electron is in the state |^(0)} initially. After the pumping 
action it is in the final state [SE'(T)), where T is the period of 
the cycle. The initial and final states |*(0)) and |*(T)} are 
related to each other through a matrix Berry phase. During 
the adiabatic cycle the electron state does not leave the doubly 
degenerate Hilbert subspace. 



A semiconductor quantum dot electron pump[l|, Q is 
a device that can make a dc current flow in unbiased 
systems. In such a device a charge can be pumped 
through the dot when the shape of the double barrier 
potential undergoes an adiabatic cyclic change. Quan- 
tum pumping is a fascinating subject and considerable 
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carried out 
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and theoretical studies have been 

ifj, El m El El, EE EE Ezl, p, 

Recently effects of noise and decoherence[21, 
2l interactions [H HE [H S H HE S3 , and 
spin[6|, |32, 123, 113] have been investigated. The pumped 
charge through a single energy level during a cycle at zero 
temperature can be calculated Q using emissivity[35[ 
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where s a p is the scattering matrix element, and a and 
j3 denote contacts 1 and 2. During an adiabatic cycle 



the parameters Ai and A2, which control the shape of the 
dot, follow a closed path C in the parameter space and 
the area defined by it is denoted as A, see Fig[T] 

The integration over the area A suggests that quantum 
pumping depends on the geometric properties of the path 
and is independent on detailed time evolution. This in- 
dicates that quantum pumping in nano-semiconductors 
may be a manifestation of a matrix (non-Abelian) Berry 
phase [IE Similar issues have been investigated in 

quite different physical systems, such as superconduct- 
ing circuits (38[. Recent investigations have shown that 
several semiconductor dots possess non-Abelian matrix 
Berry phases. Semiconductor excitons[39| and p-type 
electrons [IE Sl| can exhibit matrix Berry phases. It has 
also been shown that n-type semiconductor dots with 
spin-orbit coupling possess a matrix (non-Abelian) Berry 
phase [IE EE This is due to the presence of energy 
level degeneracy originating from time-reversal symme- 
try. 

We show in this paper that adiabatic pumping in semi- 
conductor quantum dots without spin-orbit coupling may 
be analyzed in terms of non-Abelian Berry phases. For 
non-Abelian Berry phases to be finite it is essential that 
degenerate states are present, and semiconductor quan- 
tum pumps do have non-trivial degenerate scattering 
states, see Fig[2] Each energy level of a quantum dot 
pump is doubly degenerate since there are two indepen- 
dent scattering states incident from the left and right bar- 
riers. During the pumping action these degenerate states 
can get mixed. This can be seen as follows: LetT(t+5t,t) 
be the time development operator between t + St and t, 
and let the change of the confinement potential of the dot 
during the pumping action be SV(t) = V(t + St) — V(t), 
where St is an infinitely small time interval. The degen- 
erate instantaneous eigenstates of the Hamiltonian are 
chosen as the degenerate scattering states \ipi(t)), shown 
in Figl^l Suppose that the initial electron at t is in the 
scattering state \ipi(t)). Then the overlap matrix element 
of the time development operator between the degenerate 
scattering eigenstates of the Hamiltonian is given by 



(Mt + st,t)\T(t + st,t)\Mt)) w 
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since T(t + St,t) w 1 — iStH(t). For some pumping ac- 
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tions, given by Eq. (|T9[) . this overlap is finite, and the 
state (t)) can develop into a linear combination of 
\tpi(t+5t)) and \ipi(t+St)) . This suggests that during the 
adiabatic cycle a state vector may rotate in the degener- 
ate Hilbert subspace and the initial state may not return 
to itself after the adiabatic cycle: |*(0)) ^ as 
shown in FigflJ It is unclear under what conditions this 
happens and how different the initial and final states are. 
Do they differ by a phase factor or more than a phase fac- 
tor? Theory of non-Abelian matrix Berry phases [36j is 
well suited to address these issues. 
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For a cyclic change with the period T, the instantaneous 
basis states |$i(T)) and |<I>2(T)) return to the initial 
states | $1(0)) and |$2(0)), but the coefficients ci(T) and 
c 2 (T) may not return to the initial values. In this case, 
the initial and final degenerate states may be connected 
via a finite matrix Berry phase. A matrix Berry phase is 
a non-perturbative concept and is applicable even in the 
presence of strong potential modulations that are adia- 
batic. A matrix Berry phase can arise when an energy 
level has a non-trivial degeneracy which does not split 
during an adiabatic change of the external parameters 
X p of the Hamiltonian. For a doubly degenerate energy 
level a matrix Berry phase is a 2 x 2 matrix $c connecting 
the final amplitudes (ci(T), C2(T)) to initial amplitudes 
(ci(0),c 2 (0)): 



Cl(T) )-$r( Cl{0) 
C2{T) ) c 2 (0) 



(4) 
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Using the instantaneous degenerate basis vectors the ma- 
trix Berry phase can be computed by solving the adia- 
batic version of the time-dependent Schrodinger equation 
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FIG. 2: In our model an electron moves in a one-dimensional 
asymmetric double barrier system without spin-orbit cou- 
pling. Two degenerate scattering states ■i/>i(x) and ^2(2) are 
shown. Ti, Ti and 7i,7, are wavefunction amplitudes. 



where Ay = ^^2 p (A p )i.j^- and the sum over p in Ay 
is meant to be the sum over X p . The non-Abelian gauge 
potentials are defined as 



We show that the presence of a finite matrix Berry 
phase leads to a non-zero pumped charge. Our investi- 
gation of time evolution of the electron state shows that 
the pumped charge depends on the initial state of the 
adiabatic cycle, i.e., different initial states give different 
pumped charges. Our method is applicable even when 
the potential modulations are substantial. 



(6) 



Solving the time-dependent Schrodinger equation, 
Eq.([5]), one can write the matrix Berry phase in terms 
of the non-Abelian gauge potentials (36J: 



= Pe} fcT. p A P dX P j 



(7) 



II. NON-ABELIAN U(2) GAUGE THEORY, 
MATRIX BERRY PHASE, AND PUMPED 
CHARGE 



where a path ordering P must be used since A p at differ- 
ent times are usually noncommuting. We see from this 
contour integration that 



Non-Abelian Berry phases [36| are computed from the 
exact solutions of the time-dependent Schrodinger equa- 
tion for adiabatic changes. In this approach one can 
choose the orthonormal basis states of the degenerate en- 
ergy level as the instantaneous eigenstates of the Hamil- 
tonian: H(t)\$i(t)) = E(t)\$i(t)) for i = 1,2. (Here we 
consider only adiabatic changes and other states with dif- 
ferent energies are ignored). Even when the electron state 
starts in one of the instantaneous eigenstates |$i(0)) it 
may not stay in the same eigenstate \&i(t)} at later times: 
the actual state at t may be a linear combination of the 
instantaneous eigenstates 



!*(*)) = Cl (t)|$ 1 (t)) + C2 (*)|$ 2 (t)). 



(3) 



1. Matrix Berry phase is a geometric effect, i.e. it 
depends only on the path C. 

2. Matrix Berry phase is independent of the period of 
the cycle as long as the change is adiabatic. 

It is instructive to compare this result to the Abelian 
Aharonov-Bohm phase of an electron in the electromag- 
netic vector potential A: 



e i f a A-di 



(8) 



where C is a closed path encircling the magnetic flux. 
The only difference between the Abelian and non-Abelian 



phases is the mathematical nature of the vector poten- 
tials: a number vs. matrix. Under a unitary transforma- 
tion | ip'A = J2j U*j\ipj) the non-Abelian gauge structure 
emerges 

dW 

A' k = UA k W + iU—. (9) 
For doubly degenerate levels the relevant gauge group is 

A pumped charge can be computed from the solution of 
the time-dependent Schrodinger equation, Eq.©. It can 
be calculated from the total probabilities in the regions 
to the left, right, and center of the dot, denoted by 

P h (t) = J J \*(x,t)\ 2 dx, 
P R (t) = [ L \^(x,t)\ 2 dx, 

J a/2 

Pc(t) = / \V{x,t)\ 2 dx, (10) 

J -a/2 

where L is the length of the system. The change of these 
probabilities after the nth adiabatic cycle, t — nT, is 

AP i (n)=P i (n)-P i (n-l). (11) 

for i = L,C, R. The probabilities AP^ and APr are di- 
rectly related to the pumped charge: the pumped charge 
after the first cycle is 

Q = APt(l) = -APfl(l), (12) 

provided that the change of probability in the dot, 
APc(l), is negligible. It is straightforward to relate the 
matrix Berry phase to the pumped charge using Eqs. 
(J4j) and (Tl2|) . The dependence on the non-Abelian Berry 
phases enters through Eq.JU), and the pumped charge 
depends on the initial conditions (ci(0), C2(0)) and basis 
states $j(x): 

= X;0«[(* c )«(*c)il -^^]c fc (0)*q(0), 

ijkl 

(13) 

where 

f- a/2 

O i:i = J ^(x)*^ 3 (x)dx. (14) 

III. MODEL CALCULATION OF PUMPED 
CHARGE 




FIG. 3: Transmission coefficient of a double barrier system 
with the parameters Vi/(Eoao) = 4.25 and V r /(Eoao) — 2.0. 



The simplest possible model to study quantum pump- 
ing is a one-dimensional double barrier system]^ |46|. If 
the barrier potentials are chosen to be delta functions it 
is possible to calculate analytically the eigenstates and 
the non-Abelian vector potentials. Let us consider two 
delta-function barriers with the distance a: 

V(x) = V l 6(x+~) + V r 6(x-~), (15) 

where Vi and V r are heights of the left and right barriers, 
respectively. Each energy level is doubly degenerate with 
eigenstates ipi{x) and ip2(x). in tpi(x) (tp2(x)) the elec- 
tron is incident from the left (right), see Fig.l. Due to 
time reversal invariance the scattering amplitudes satisfy 
T3/V1 = t^/t[. In the following we consider the pumped 
charge of an electron in the energy level E — where 

is a wavevector and m is the electron effective mass. 

r 2 

We use the energy unit Eq — 2 ^ ai and the length unit 
ao = 2a/3. 

In the numerical calculation of the pumped charge it 
is convenient to use box normalization with length L. 
Let us comment on how the limit L — > 00 is taken. It 
is natural to choose the incoming scattering states from 
the left and right as the instantaneous basis states: ipi(x) 
and %l>2{x) (any other orthonormal set can also be used 
to calculate the non-Abelian vector potentials). However, 
for a finite value of L these two states are not orthogonal. 
We choose one of the instantaneous basis state as tpi(x) 
and construct the other basis state via Gram-Schmidt 
orthogonalization using ip2 (x) : 

$i(x) = ipx{x), 

$ 2 (a:) =aMx)+0Mx)> ( 16 ) 

As L — * 00 the overlap (VulV^) — * and $2(2:) — * ifoix). 
Below the pumped charge is calculated for sufficiently 
large values L where the size dependence is negligible. 
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A. Parity symmetry: absence of matrix Berry 
phase and pumped charge 

When the double barrier potential retains parity sym- 
metry during the adiabatic cycle, i.e. Vi(t) — V r (t), the 
pumped charge is zero. (The adiabatic parameters are 
(Vi(t),a(t))). This follows from the expression for the 
pumped charge given in Eq.([T]): we have verified that 

the imaginary part of ^r^ 1 + is zero. We 

show that this is fully consistent with the absence of a 
matrix Berry phase. For this purpose it is convenient 
to choose the following orthonormal instantaneous basis 
states: 



*i(a:) - 
$2 (a) = 



1 



/]VT 
i 

m 



:(lpl(x) ~ lfo{x)), 



(17) 



where the normalization factors are N\ = 2 + 2 (V'l I ^2) 
and N2 = 2 — 2(ipi\ip2)- We will use this orthonormal set 
to calculate the non-Abelian vector potentials [47]. Since 
two barrier heights are equal, i.e. Vi = V r , it can be 
easily shown that ipi( x ) = V J 2( — x). From this it fol- 
lows that $i(x) and $2(2;) are, respectively, even and 
odd functions in x. Taking the above states as basis vec- 
tors, and using the results (V>i| gfHV'i) = W^lgf-I^a) 



and 



i^Apiv^/ — \vi\ alp 1 vi/ w it n the orthogonality 
($i|$2) = we find that the following off-diagonal ele- 
ments of vector potentials are zero: 



((Vil^-l^ + ^il^-l^) 



dXr, 



Similarly {A p )2i = can be proved. Since the off- 
diagonal elements of the vector potentials are zero the 
matrix Berry phase is zero. This effect of parity sym- 
metry is already known from the investigation of matrix 
Berry phases of isolated II- VI and III-V quantum dots 
in the presence of spin-orbit coupling: when the in-plane 
confinement potential of the dot has parity symmetry the 
matrix Berry phase vanishes [43. l43l|. 



B. Numerical calculation of pumped charge 

When the parity symmetry is broken during the adi- 
abatic cycle, i.e. Vi(t) 7^ V r (t), a pumped charge is ex- 
pected. We choose the adiabatic parameters as Vi and 
V r and select elliptic paths 



{V l (t),V r {t)) = 

(Vi, c + AV/ cos(wt), K,, 
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AV r sm(u)t)). (19) 
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FIG. 4: Numerical results for the change in the total probabil- 
ities in different space regions after the first adiabatic cycle as 
a function of wavevector: APl(1) (squares) and APr(1) (cir- 
cles). The pumped charge Q/e calculated analytically using 
Eq.([lJ (triangles). We have used L/ao = 500 and w/Eq = 0.1. 



At t — we start with the scattering state \ipi). Wc 
investigate pumped charges when the potential modula- 
tions are substantial : AVi/Vi lC — 0.7, AV r /V 7%c = 0.7, 
where Vi >c /(a E ) = 2.5 and V r;C /(a E Q ) = 2.0. At the 
start of the adiabatic cycle the quantum dot is charac- 
terized by the scattering amplitudes shown in Fig[3] The 
time dependent Schrodinger equation, Eq.© is solved 
numerically, and we determine the amplitudes Ci %(T), 
from which the matrix Berry phase can be computed 
using Eq.((4]). The pumped charge after the first cycle, 
APl(I), is calculated as a function of wavevector using 
Eqs.|in]). The results are shown in FigQJ The pumped 
charge is maximum near kao = 1.25, which is close to 
the position of the first resonance in the transmission 
coefficient, ka = 1.56[45|, [46J. Note that the vanished 
probability on the right side of the dot appears on the 
left side: AP L (1) = -AP R (1). This implies that the 
pumped charge does not accumulate in the dot, which is 
consistent with Eq.|T|), see Ref.Q- We have also com- 
pared quantitatively our results for the pumped charges 
with those obtained from Eq. (JTJ) : they are in a good qual- 
itative agreement. 

In our approach we can follow microscopically how 
various quantities actually depend on time. Fig[5] dis- 
plays time dependence of the total probabilities in dif- 
ferent space regions during consecutive adiabatic trans- 
formations. Results are shown for kao = 1-30. We note 
that Pl(ji) for n = 0, 1, ... are all different, implying that 
pumped charge APl(u) depends on the number of con- 
secutive adiabatic cycles. This is because initial state of 
the nth cycle fy((n — 1)T) is different from that of the 
(n - l)th cycle *((n - 2)T). In Eq.fT]) this dependence 
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t/T 



FIG. 5: Time dependence of the total probabilities in different 
space regions, Pl (solid) and Pr (dashed), during consecutive 
adiabatic transformations. 



final states are connected via a finite matrix Berry phase. 
This difference between initial and final states of an adia- 
batic cycle implies that a pumped charge is present. The 
main physics behind this is that even when the electron 
state starts in one of the instantaneous eigenstates it may 
not stay in the same eigenstate at later times: often the 
actual state at time t turns out to be a linear combination 
of the instantaneous eigenstates. One can prepare the ini- 
tial state to be the incident scattering state from the left. 
When consecutive adiabatic cycles are performed start- 
ing from this initial state the pumped charge at the nth 
cycle is not necessarily the same as that of the (n — l)th 
cycle. Our investigation of the time dependence thus 
shows that the pumped charge depends on the nature of 
the initial state of a cycle. Moreover, during the period 
of a cycle the probability flux through the dot displays a 
significant variation in time. It would be interesting to 
test these results experimentally. 
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